First-principles modeling of electrostatically doped perovskite systems 
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Macroscopically, confined electron gases at polar oxide interfaces are rationalized within the sim- 
ple "polar catastrophe" model. At the microscopic level, however, many other effects such as electric 
fields, structural distortions and quantum-mechanical interactions enter into play. Here we show how 
to bridge the gap between these two length scales, by combining the accuracy of first-principles meth- 
ods with the conceptual simplicity of model Hamiltonian approaches. To demonstrate our strategy, 
we address the equilibrium distribution of the compensating free carriers at polar LaAlOs/SrTiOa 
interfaces. Remarkably, a model including only calculated bulk properties of SrTiOs and no ad- 
justable parameters accurately reproduces our full first-principles results. Our strategy provides a 
unified description of charge compensation mechanisms in SrTiOs-based systems. 

PACS numbers: 71.15.-m, 71.70.-d, 73.20.-r 
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The unusual behavior of the (001) LaAlOa/SrTiOa 
interface is commonly understood in terms of the "po- 
lar catastrophe" (PC) model. [1] By stacking charged 
(LaO)+ and (A102)^ layers on top of the charge- neutral 
Ti02 and SrO layers, one obtains a net interface charge 
density of apc = +e/2S, where S is the unit-cell cross 
section. This produces a diverging electrostatic energy, 
unless apc is neutralized by an external free charge, 
which would explain the appearence of confined mobile 
carriers at this interface. 

This model, while appealing, misses many important 
effects, that are crucial for a realistic description of the 
interface. For example, it was shown that strong polar 
distortions in LaAlOa Q and/or in SrTiOa 0,01 partially 
screen the excess charge, delaying the onset of metallic- 
ity far beyond what the PC arguments would predict. 
Next, it was shown that H adsorbates [HI or oxygen va- 
cancies at the open LAO surface can profoundly alter 
the sheet density of free carriers. Moreover, reversible 
metal-insulator transitions can also be induced upon ap- 
plication of an external bias 0, § . Both effects go clearly 
beyond the oversimplified PC description. Finally, the 
PC model cannot predict truly microscopic properties 
of the system, such as the spatial decay and confine- 
ment of the free electrons near the interface. In an at- 
tempt to answer these latter crucial questions, various 
quantum-mechanical explanations were proposed , 
but their relative importance, especially in relationship 
to the macroscopic electrostatics arguments, is unclear. 

A route towards bridging the gap between classical 
electrostatics and quantum theory in LAO/STO and re- 
lated systems was proposed in Ref. S The strategy is 
based on the "formal" |lll| definition of the polarization, 
P, in quantum-mechanical systems, which has a simple 
classical interpretation in terms of a point-charge model 
(see Fig. [1]). The dipole moment of an individual (LaO)- 
(AIO2) unit is d = —ea/2 [black arrows in Fig. IHb)], 
where a is the out-of-plane lattice parameter and e is 
the (positive) electronic charge. This corresponds to a 
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FIG. 1. a) Scheme of the Ti02:LaO interface, showing the 
built-in (P", black arrows) and induced (AP, blue arrows) 
polarization, and the free charge crfi-cc- b) Electric displace- 
ment as a function of z, illustrating Eq. [T] c) Calculated 
internal electric fields £ in bulk STO and LAO as a function 
of D. The polar nature of the interface stems from the impos- 
sibility of finding a common value of D for which ^lao and 
SsTO are both zero. 

"built-in" polarization Plao ~ "6/25, where S is the 
cell surface. Conversely, -Pgxo ~ ^ because the STO lay- 
ers are formally charge-neutral. Note that there is no 
left-over ionic charge at the interface - we have reinter- 
preted (Tpc as a surface density of bound charge, that 
arises because of a discontinuity in P. 

Depending on the electrical boundary conditions, 
macroscopic electric fields, respectively £lao and fsTOi 
will be present on one or either side of the junction. The 
electric fields will in turn perturb the individual LAO and 
STO layers, producing an "induced" polarization that we 
call APlao and APgxo ■ If we now define the total po- 
larization as P = P°-t- AP, and the electric displacement 
as D = eoS + P, an exact relationship follows, Q 
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where CTftoc is a surface density of "free" charge con- 
fined to the interface region. Eq. [T] generalizes the PC 
model by taking rigorously into account the effect of po- 
lar distortions (AP is implicitly contained in D), ex- 
ternal biases [£{D) is a bulk property of either mate- 
rial, and is a unique function of D, see Fig. [UJc)] and 
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charged species adsorbed on the far-away surfaces (the 
flux of D corresponds to the surface charge density). By 
appropriately choosing the two independent parameters 
-Dlao and Dsto one can therefore describe the local 
properties of an ideal interface within arbitrary bound- 
ary conditions, encompassing virtually all theoretical ap- 
proaches (stoichiometric or non-stochiometric superlat- 
tices and various flavors of slab geometries) that were 
used so far in the literature. 12, ll3| One does not need 
to worry about the specific mechanisms and/or supercell 
geometries that determine a certain equilibrium value of 
-Dlao and Dsto, as long as the interface can be thought 
as isolated (say, separated by at least two or three unit 
cells of LAO and STO on either side) . 

To work our way towards the microscopies, it is now 
tempting to take the analogy to macroscopic Maxwell 
equations one step further, and write 



dD{z) 
dz 



Pfrco(z). 



(2) 



Here phcc{z) is the spatially resolved planar average of 
free carriers, whose integral along z yields CTfroo- It is 
easy to verify that Eq. [5] is consistent with Eq. [T] Here 
one runs into trouble, however, as one needs to estab- 
lish a truly microscopic definition of both D{z) and pbcc- 
(Note that this is not necessary at the level of Eq. [l] 
which deals only with macroscopic quantities.) This is a 
nontrivial issue in a typical metal, where the polarization 
(and hence D) is ill-defined. Furthermore, phee is micro- 
scopically difficult to identify, as the bands correspond- 
ing to the conduction electrons are generally entangled 
with lower-lying bound states. In a doped oxide or semi- 
conductor, however, the valence and conduction bands 
usually preserve their identity, i.e. a well-defined energy 
gap persists between conduction-band and valence-band 
states. This naturally leads to a definition of pfree based 
on the overall density of the partially occupied states 
near the Fermi level. The remainder is an integer num- 
ber of electrons that we identify as bound charges. We 
use these latter orbitals to define a layer-resolved electric 
displacement based on a Wannier decomposition of the 
polarization, [3] in analogy with standard insulators. 

We are now ready to verify Eq. [2] directly on our first- 
principles calculations. To provide a representative num- 
ber of test cases, we study three combinations of J5sto 
and I?LAOj which are summarized in the insets of Fig- 
ure [21 Case (a) corresponds to full compensation, e.g. 
that of a thick LAO overlayer on a thick STO sub- 
strate. Case (b) corresponds to partial compensation, 
which can occur at intermediate LAO thicknesses [l5| . 
or in the case of an electrical bias applied between the 
electron gas and an electrode deposited on the free sur- 
face. Case (c) physically corresponds to a "back-gating" 
regime, where an electrical bias is applied between the 
electron gas and an electrode placed at the other end 
of the STO substrate. In practice, we use slab geome- 
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FIG. 2. Conduction charge (a-c) and local electric displace- 
ment (d) at the SrTiOs/LaAlOs interface. In (a-c) the solid 
curves are phcc{z) and the dashed curves are the nanos- 
moothed pfrce(z). Insets show an approximate diagram of the 
effective local potential. In (d) the symbols show the Wannier- 
based local electric displacement computed from the bound 
charges. The curves represent phee{t)dt + Dsto- 



tries of the type vacuum/ (SrTi03)„/(LaA103)m/vacuuni 
(we use n — IQ and m = 3 in our calculations), where 
the boundary conditions on D are enforced as explained 
in the Supplementary Information. (All the other rele- 
vant computational parameters are also described there.) 
In Fig. [2] we show the relaxed phoc{z) (a-c) and the 
layer-by-layer (locally averaged) electric displacement Di 
(d) for each combination. In Fig. [D^d) we also plot 
three curves that we constructed by numerically inte- 
grating the nanosmoothed charge densities, pi-cee{z)- [We 
are therefore verifying the integral version of Eq. [2l 
D{z) = f^^Phee{t)dt + D(— oo).] The matching is ex- 
cellent in all cases, demonstrating the high accuracy of 
Eq. [21 Note that in one of the examples (case c) pfree de- 
cays to zero relatively fast when moving away from the 
interface, while it spreads over the whole volume of the 
SrTiOa film in the other two cases. This is due to the fact 
that in c) the asymptotic electric field in STO is not zero, 
but equal to £{D = -0.2) —12 MV/m [see Fig. [flc)]. 
This produces a confining wedge potential, that limits 
the spread of pfico- Conversely, in (a) and (b), £ vanishes 
at z — — oo, and the outermost electrons are only loosely 
bound. 

Eq.[2]is an important result, in that it establishes a di- 
rect, virtually exact relationship between the density of 
compensating carriers and the local polarization in LAO 
and STO. This answers pressing experimental questions 
concerning precisely this point, as_ polar distortions in 
SrTiOa were recently observed. 
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found implications over the theoretical understanding of 
electron confinement in this system, as we shall demon- 
strate in the following. 

Essentially, the equilibrium distribution of the con- 
duction charge is determined by two competing effects. 
One is the electrostatic energy, that tends to localize 
the electrons as close to the interface as possible. The 
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FIG. 3. Bulk properties of SrTiOs used in the tight-binding 
model. Upper panels: internal potential (a) and dielectric 
constant (b) as a function of D. Lower panels: t2g conduction 
band structure for D = (c) and D = —e/2S (d). 

strength of the attraction will depend on the static di- 
electric constant of the underlying insulator. The other 
is the quantum-mechanical kinetic energy of the elec- 
trons. This will tend to spread the electrons in space, 
with a strength that depends on the band dispersion. To 
see whether, and to what extent, the large polar distor- 
tions in STO affect these competing driving forces, we 
performed calculations of bulk SrTiOa as a function of 
the electric displacement [Tt'I, by covering the range of 
^STO values that are relevant for the LAO/STO system. 
For each value of D we extract the built-in electric field, 
the total internal energy and the relevant parameters of 
the lower part of the conduction band. These are the 
tight-binding hopping integrals between Ti-derived or- 
bitals with t2g symmetry {dxy^ d^z and dyz)- As the t2g 
orbitals are fairly well localized in space, it is sufficient 
for the present study to consider only the first three shells 
of nearest-neighbors Ti sites. 

In Fig. inja) we show the electric field as a function of 
the displacement field D. Note the strong nonlinearity, 
which is evident in the plot of the dielectric constant as a 
function of D [Fig. [Hb)]. In Fig. [JJc) we show the band 
structure as it results from the third-neighbor Hamilto- 
nian, for the centrosymmetric cubic state at _D = 0. Note 
the symmetry of the bands, which are characterized by a 
three-fold degeneracy at F. A polarization [the extreme 
case D = —e/2S is shown in Fig.jSjd)] lifts this degener- 
acy, by producing a strong splitting at F between the de- 
generate dxz/dyz orbitals and the d^y orbital. This split- 
ting is dominated by the strong reduction in the dxz/dyz 
bandwidth along the T ^ X and T Z directions - the 
corresponding hopping terms are reduced by as much as 
30% and 25%, respectively. Polarization-related changes 
in other matrix elements appear to be less pronounced. 

We shall now use these data to develop a quantitative 
model of the equilibrium distribution pt].ee(-z)- We make 
a rather bold assumption here, and state that the role 
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FIG. 4. Electron confinement and spatial distribution. Red 
curves with square symbols are the first-principles results, 
black curves (circles) are the results of the model. Param- 
eters are the same as in Fig. [2] (a-c). (d) shows a thicker 
24-ceU STO film. The dashed (blue) and dot-dashed (green) 
curves were obtained by artificially altering selected features 
of the model (see text). Insets: tight-binding band structures 
corresponding to the green (di) and black (d2) curves in (d). 
The arrow indicates the lowest dy^/d^x state. 

played by the LaAlOa overlayer in determining p{^ee{z) 
is marginal, except for two crucial effects: i) it con- 
fines the conduction electrons to the STO side, and ii) 
it defines the electrical boundary conditions through the 
value of Dlao- Based on this Ansatz, we represent the 
LAO/STO interface systems discussed in the previous 
paragraphs as pure STO slabs, periodic in plane and n- 
layer thick, where the boundary values of the electric 
displacement field at the two surfaces are set to -Dsto 
and Z^LAO- To each Ti site I we assign three orbitals of 
t2g symmetry, and a charge density pi. The charge den- 
sity defines the local value of the electric displacement Di 
through Eq. [2l The Hamiltonian matrix elements are de- 
fined by the electrostatic potential Vi [calculated from Di 
using the bulk Vsto(-D) of Fig. [2]Ja)], which rigidly shifts 
the on-site terms, and by the Dj-dependent hopping pa- 
rameters that we interpolate from the bulk SrTiOs data. 
Upon diagonalization we obtain the wavefunctions, that 
self-consistently determine pi within the constraint Eq. [T] 
In Fig. m we compare the results of the model to the 
first-principles simulations discussed earlier. We include 
in the comparison a fourth simulation that we did for 
case I?sTO = 0, I^lao = — e/25, but with a thicker STO 
layer [n — 24). The agreement is remarkably good. This 
indicates that the polarization-dependent bulk proper- 
ties of SrTiOa, together with the boundary values of D, 
are sufficient to explain the distribution of conduction 
charge in this system. This suggests that the interaction 
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between Ti- and La-derived orbitals is not an essential 
factor in determining electron confinement, contrary to 
the conclusions of Ref. 0. Binding of the electrons to the 
interface is indeed guaranteed by Eq. [21 

Now that we have a reliable model we can directly 
quantify the impact of each specific STO bulk property 
on the distribution of phoc- First, if we neglect the non- 
linearity in the dielectric permittivity esTo(-D)) and in- 
stead use a constant esTo(-D) — esTo{0) ~ 500 we obtain 
a much broader distribution [blue curve in Fig. El^d)]. 
This indicates that the carrier distribution is strongly 
sensitive to the dielectric properties of bulk STO; this 
seems to be an accepted fact in the experimental com- 
munity, [3| but has received surprisingly little attention 
in earlier ab-initio studies. Second, if we suppress the 
D-dependence of the STO band structure, and use the 
D = t2g Hamiltonian throughout the film, we obtain 
[green curve in Fig. IH^d)] an excessive accummulation of 
charge in the near-interface region. This effect can be 
understood by comparing the self-consistent bandstruc- 
tures of the original [Fig. |4i;d2)] and the "i2g(-D = 0)" 
[Fig.mjdi)] tight-binding models. In both cases there is a 
strong splitting at F between the dxy and d ^^/dy ^ bands, 
in agreement with the findings of Refs. [ij, [ij and [l^- 
In Fig. Illjdi), however, this splitting is only induced by 
confinement effects due to the wedge-like electrostatic po- 
tential near the interface po| . The polarization-induced 
perturbations in the STO t2g bands [Fig. [3jd)] signif- 
icantly enhance such a splitting [Fig. |4l^d2)] and shift 
the dxz/dyz bands further up in energy. (The effect is 
strongest on the lowest dxz/dyz band, marked with a 
red arrow in the figure). This upshift, in turn, pushes 
the weight of the d^z/dyz electrons away from the LAO 
interface, which explains the difference between the re- 
spective electron distributions [green and black curves 
in Fig. lU^d)]. Note that an analogous t2n splitting was 
experimentally observed in LAO/ STO, [21] and theoret- 
ically also discussed in the context of the closely related 
LaTiOg/SrTiOa system. [H 

The tight-binding method used here has clear points of 
contact with the strategy of Refs. 22 and 23,. However, in 
our approach there is a crucial innovation. Here, at dif- 
ference with Ref. [22, we extract all the ingredients of the 
model from bulk calculation of pure insulating SrTiOa, 
without including any adjustable parameter. This forces 
us to build a universal and transferable model, which can 
be readily applied to essentially any situation involving 
electrostatic doping of SrTiOs, and is not restricted to 
the specifics of the LaAlOs/SrTiOa interface. For exam- 
ple, our strategy could be readily used, with little modifi- 
cations, to interpret the recent findings of electron gases 
at the bare SrTiOa surface. [20] More importantly, our 
model could be readily extended to account for other 
physical ingredients not considered here, e.g. strong cor- 



we take as input. This is an enormous advantage, both 
conceptually (the model is based on few parameters that 
are easy to interpret) and practically (the tight-binding 
model is several orders of magnitude more efficient than 
a full first-principles calculation). More generally, our re- 
sults open exciting new avenues for the study of confined 
electron gases in oxide systems, with optimal accuracy 
and dramatically reduced computational cost. 
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relations 'l2'| and strain effects; 2J| all we need to do 



is to refine the theoretical description of bulk STO that 
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